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l. Übersicht 
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Um eine Aussage über das Verhalten eines geplanten Windrades zu er- 

halten oder um Grundlagenforschung an Windrädern zu betreiben, gibt 

es drei prinzipiell unterschiedliche Möglichkeiten: 

a) Man versucht, auf rein theoretischem Wege Ergebnisse zu 
entwickeln. /15/. 

b) Man führt Versuche im Windkanal durch. 

c) Man untersucht bereits gebaute Windräder und versucht die hier 
gemachten Erfahrungen auf das eigene Problem 
zu übertragen. 

Mit der Problemstellung b) und c) setzt sich diese Studienarbeit 


auseinander. 


Da im Windkanal ganz andere Strömungsverhältnisse vorliegen als 
"draußen", muß man sich zunächst Gedanken darüber machen, wie das 
getestete Windrad sich außerhalb des Kanals verhalten würde. Dies 
geschieht in Abschnitt 7. Liegen diese Ergebnisse vor oder unter- 
sucht man, wie unter c), ein bereits gebautes Windrad, so können 
die Daten auf ein Windrad anderer Größe umgerechnet werden. Eine 
andere Möglichkeit besteht darin, die Umrechnung auf andere 
Strömungsverhältnisse auszurichten. Diese Rechnungen werden in 
Abschnitt 8 durchgeführt. Dabei geht es vor allem um dimensionslose 
Kenngrößen, insbesondere um Wirkungsgrade. Daß sich nicht nur die 
Wirkungsgrade bei einer anderen Windradgröße verändern, sondern 
auch bei der Berechnung der Gesamtwindleistung einiges zu beachten 
ist, wird in Abschnitt 9 gezeigt. Doch zunächst zu den einleitenden 
Abschnitten... 


2. Einleitung 
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Bevor der Mensch lernte, sich seine Umwelt zunutze zu machen, war 
er nur auf seine Muskelkraft angewiesen. Dann begann er das Feuer 


zum Heizen und Tiere für verschiedene Arbeiten einzusetzen. 


Ab wann der Mensch den Wind für sich arbeiten ließ, ist unbekannt. 
Die Ägypter z.B. benutzten das Segel schon 2800 v.Chr. Vor mehr als 
3000 Jahren drehten sich dann schon die ersten Windräder. Im euro- 
päischen Raum verbreitete sich zunächst eine Bauart mit einem stei- 


nernen Turm, auf dessen Spitze sich ein mehrflügliges 


Windrad drehte. Das Windrad selbst bestand aus einem System von 
"Dreieckssegeln". Später baute man Windräder mit vier Flügeln, die 
aus einem stoffbespannten Holzgerippe aufgebaut wurden. 


Die erste entscheidende Verbesserung brachte das Spätmittelalter in 
Gestalt der "Deutschen Bockwindmühle". Bei diesem Typ konnte erst- 
mals das auf einem Bock drehbar gelagerte Mühlenhaus mit dem Wind- 
rad in die richtige Windrichtung gedreht werden (Abb. 1). 


In verschiedenen Schritten wurde die sich rasch verbreitende Bock- 
windmühle vor allem in Holland verbessert. Mit der sog. "Holländer- 
windmühle", bei der nur noch die Turmkappe mit dem Windrad gedreht 
wurde, während das Mühlenhaus feststand, erreichte die historische 
Windmühle ihre Vollendung (Abb. 2). 


In den Vereinigten Staaten wurde im 19. Jahrhundert die "Amerika- 
nische Windturbine" (Abb.3) auf entlegenen Farmen in großer Zahl 
eingesetzt. Zunächst wurde sie zum Wasserpumpen gebraucht, später 
aber auch zur Versorgung der Farmen mit elektrischer Energie. 


Seit man begann das Feuer auch zur mechanischen Energiegewinnung 
mittels Dampfmaschinen zu nutzen, verringerte sich die Bedeutung 
der Windkraft mehr und mehr. Zur Stromerzeugung ‚schien die Wind- 
kraft ungeeignet und unwirtschaftlich zu sein, denn es hatten sich 
große Kohle-, Öl- und Atomkraftwerke zusammen mit dichten Versor- 
gungsnetzen entwickelt. 


Erst in neuerer Zeit kam die Windenergie wegen der Suche nach rege- 
nerativen und umweltfreundlichen Energieformen wieder in das 
öffentliche Interesse. Dabei kann man heute auf Erfahrungen zu- 
rückgreifen, die auch in Zeiten allgemeinen Desinteresses an der 
Windenergie gemacht wurden. Die Windmühlen von damals haben sich 
dabei zu Windturbinen gewandelt, die - wegen der Anwendung von mo- 
dernen Aerodynamikkenntnissen - nur noch wenig Ähnlichkeit mit 
ihren Vorfahren haben. 
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Abb, 1 


Deutsche Bockwindmühle 











Abb. 2 


Holländerwindmühle 











Abb. 3 


Amerikanische Windturbine 








3. Einteilung von Windkraftanlagen 
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Windkraftanlagen werden in der Hauptsache nach folgenden Merkmalen 
eingeteilt: | 


3.1 Arbeitsprinzip 


a) Widerstandsprinzip: Der Wind treibt eine luftundurchlässige 
Fläche vor sich her. Die treibende Kraft ist die Widerstandskraft. 


b) Auftriebsprinzip: Die Energie wird der Luftströmung durch einen 
mehr oder weniger aerodynamisch gestalteten Körper entzogen, der so 
zum Wind angestellt ist, daß die Auftriebskraft die Anlage an- 
treibt. Windräder, die nach dem Auftriebsprinzip arbeiten, haben 
einen höheren Wirkungsgrad. 


3.2 Achslage: 


Man unterscheidet: 
Horizontalachsrotor, Wagner-Rotor und Vertikalachsrotor. Vertikal- 
achsrotoren haben den Vorteil, daß sie nicht zur Windrichtung 


eingestellt werden müssen. 


3.3 Schnellaufzahl (siehe 5.2) 


Schnelläufer: Schnellaufzahl > 3 
Langsamläufer: Schnellaufzahl X 3 
Schnelläufer erzielen höhere Wirkungsgrade. Ihr Anlaufverhalten ist 


schlechter. 
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3.4 Blattzahl 


Es gibt den 1-, 2-, 3-Blattrotor ... und die Windturbine. Eine hö- 
here Blattzahl verursacht eine Verminderung der Schnellaufzahl. 


In jüngerer Zeit scheinen folgende Begriffe zur Beschreibung der 
Größe einer Windkraftanlage gebräuchlich zu werden: GROWIANE (Groß- 
windanlagen), Leistung im MW-Bereich, Rotordurchmesser um 100m; 
KLEWIANE (Kleinwindanlagen), Leistung im KW-Bereich; MIWIANE (Mitt- 


lere Windanlagen). 


Weitere Merkmale ergeben sich aus besonderen konstruktiven Ausfüh- 


rungen. 


3.5 Beispiele für Windräder 


Leistungsbeiwert c_ (siehe 5.2) ca: 


P 
Schalenkreuz - 
Savonius-Rotor 0,2 
Darrieus-Rotor 0,38 
Zweiblatthorizontalachs-Rotor 0,48 
Windturbine 0,36 
Sonderbauformen: 


Savonius-Darrieus-Rotor 
Mantelturbine 
Doppelrotor 
Konzentrator 


Skizzen dieser Windräder sind in Abb. 4 gegeben. 
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2 Blatt-Rotor cyp=0,48 Kl 
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Amerikanische wWindturbine 
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Darrisus-Rotor 
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Savonius-Darrisus-Rotor 


I 


wWindrad mit Konzentrator 





en =0;,48 
cp’=1,15 (A, -RotBrfläche) 


Abb. 4 Beispiels für Windräder 


4. Der Wagner-Rotor 


Durch den Wagner-Rotor wird die Idee aufgegriffen, den Achswinkel 
deutlich verschieden von 0° bzw. den Konuswinkel deutlich verschie- 
den von 90° zu wählen (Abb. 5). Dabei ist zumindest die schräge 
Achse schon ein Stück Geschichte: Wie Abb. 6 zeigt, hatte der soge- 
nannte Tjasker aus dem 16. Jahrhundert bereits eine schräge Achse. 
Die Verlängerung der Achse ist eine "Archimedische Schraube", die 
als Wasserpumpe dient. 


Da die Energiedichte des Windes gering ist, wurden verschiedene 
Versuche unternommen, um diese zu erhöhen. Das kann z.B. durch eine 
Ummantelung eines Windrades geschehen (Abb. 4) oder durch eine ge- 
zielte Ausnutzung von Wirbelbildungen /5/. Bei der Suche nach gün- 
stigen Standorten für die Windenergienutzung wird darum auch der 
Offshore-Bereich mit berücksichtigt, weil hier schon bei niedrigen 
Höhen hohe Windgeschwindigkeiten vorherrschen. Darüber hinaus weht 
der Wind über See wesentlich gleichmäßiger als über dem Land. Will 
man nun die Windenergie auf See nutzen, so ist man mit Horizontal- 
achsrotoren auf einen Turm angewiesen, der, wie es bei Leuchttürmen 
vor der Küste der Fall ist, am Meeresboden fest verankert sein muß. 
Eine Montage solcher Großwindanlagen auf einem Schiff scheint nicht 
möglich, weil das Schiff durch die Kräfte am relativ hohen Turm 
kentern würde. Durch den Wagner-Rotor besteht nun die Möglichkeit, 
wie in Abb. 7 gezeigt, auch eine Großwindanlage auf einem Schiff zu 
montieren und die hohe Energiedichte auf See auszunutzen, weil man 


bei geschickter Winkelwahl keinen Turm benötigt. 


Der Wagner-Rotor arbeitet nach dem Auftriebsprinzip. Es handelt 
sich um einen Schnelläufer. Blattzahlen von 1 1/2, 2 und 4 werden 
angewandt. Der Wagner-Rotor muß der Windrichtung nachgeführt wer- 
den; dabei steht er so zum Wind, daß die Luft in den durch die Flü- 
gel gebildeten "Trichter" hineinströmt. Der Wagner-Rotor eignet 
sich sowohl als Großwindanlage als auch als "KLEWIAN". 
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Per Wagner - Rotor 


Konuswinkel 


Er IN Achswinkel 


Abb. 5 Winkel an asner-Rotor 
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Abb. 7 Wagner-Rotor, Land und See /11/. 


—+- - 
- a 





- n 2 7 le on 5, - 
= ER * “ + Ba . es x = Ge 
a er en Te 





nen N 


ie mine ya 


T 
9’ 
! 


TEE NER RR N 


5. Gleichungen zur Windradberechnung 
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5.l Berechnung von Windrädern 


Die Formeln zur Windenergie werden aus den Grundgleichungen der 
Strömungsmechanik gewonnen: Kontinuitätsgleichung, Energie- 
erhaltungssatz und Impulssatz. Es handelt sich zunächst um eine 
eindimensionale Betrachtung einer inkompressiblen, reibungsfreien 
Strömung. Mit Hilfe der empirisch ermittelten Daten Gleitzahl £ und 
Auftriebsbeiwert Ce, können die wirklichen Verhältnisse am Windrad 
berücksichtigt werden. Durch die Gleitzahl und den Auftriebsbeiwert 
geht eine Aussage über Zähigkeit und Turbulenzgrad in die Gleichun- 
gen der Wirkungsgrade ein, die das Windrad beschreiben. Die indu- 
zierte Wirbelbildung, die vom Druckunterschied an der Profilober- 
und -unterseite an den Flügelenden herrührt, wird durch einen ge- 
sonderten Wirkungsgrad, den Blattzahlverlustfaktor 7z berücksich- 
tigt. Die Eindimensionalität wird durch eine Integration über die 
ganze Windradfläche überwunden. Man beachte jedoch, daß die Wind- 
radberechnung weiterhin eine inkompressible Strömung voraussetzt. 
Genaueres zur Windradberechnung kann dem Seminarvortrag, /15/, ent- 
nomnen werden. 


Da es in diesem Rahmen nicht möglich ist auf die grundlegenden For- 
meln zur Windradberechnung genauer, etwa mit Beweis, einzugehen, 
soll hier nur eine Auflistung der wichtigsten Formeln gegeben wer- 
den. Man kann folgende Literaturquellen für ein genaueres Studium 
verwenden: /6/, /7/, /8/, /9/, /ı0/, /11/. Formeln, die hier nicht 
zu finden sind, werden im Anhang hergeleitet. 

In Abb. 8 ist die Strömung an einem Flügelelement mit den dort auf- 
tretenden Geschwindigkeitskomponenten gezeigt. Abb. 9 erklärt die 
Winkel am Flügel: Windangriffswinkel, Blattwinkel und Windeinfalls- 
winkel. Die wichtigen Einzelwirkungsgrade am Windrad sind in ihren 
Abhängigkeiten in Abb. 10 dargestellt. 


5.2 Definitionen 
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5.3 Formelsammlung 
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8 Strömung am Flügelelement 
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Abb. 9 Winkel am Flügel /11/ 
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Abb. 10 Darstellung der Wirkunzsgrads 
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6. Strömungsmodelle und Ähnlichkeitstheorie /6/ 
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Das Studium dynamisch ähnlicher Strömungen bildet die Grundlage für 
Modelluntersuchungen. Strömungsvorgänge um Flugzeuge, Schiffe, in 
Flüssen und nicht zuletzt auch um Windräder sind recht komplex, so 
daß oft rechnerisch nur überschlägige Aussagen möglich sind. Durch 
Untersuchungen der geplanten Großausführung in Modellen versucht 
man verbesserte Aussagen zu erhalten. 


Geometrische Ähnlichkeit beinhaltet das Verhältnis korrespondieren- 
der Längen in Natur und Modell. Der Längenmaßstab ist: 
Im 
I, Gl. 30 
N 

Dynamische Ähnlichkeit ist vorhanden, wenn Natur und Modell bei 
vorgegebener geometrischer Ähnlichkeit zu entsprechenden Zeiten 
ähnliche Strömungsverhältnisse aufweisen. Dies ist der Fall, wenn 
das Übertragungsverhältnis der Kräfte, die auf korrespondierende 
Fluidteilchen in der Natur und im Modell wirken, für jede Kräfteart 
die gleiche ist. Dabei kommen insbesondere die folgenden Kräfte- 
arten in Betracht: 

Trägheitskräfte 

Schwerekräfte 

Zähigkeitskräfte 

Druckkräfte. 
Ein konstantes Übertragungsverhältnis der einzelnen Kräftearten 
bedeutet, daß auch die Quotienten zwischen zwei verschiedenen Kräf- 
tearten in Modell und Natur gleich sind. Zu berücksichtigen sind 
daher Maßstäbe für Längen, Zeiten, Kräfte, Geschwindigkeiten, Be- 
schleunigungen, Massen, Volumina, Flächen, Dichten usw. Jedoch kön- 
nen nur drei von ihnen unabhängig gewählt werden. Wir führen daher 
neben dem Längenmaßstab den Zähigkeitsmaßstab und den 
Geschwindigkeitsmaßstab ein: 


A Gl... 31 


„_M 
Ay v. Gl. 32 
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Sowohl in der Natur als auch im Modell müssen die Navier-Stokes- 
Gleichungen erfüllt sein. Macht man diese dimensionslos, so erhält 


man die Größen 












Rt „ Irägheitskräfte „ ReyNOLDS-Zahl ° Gl. 33 
oe v Zähigkeitskräfte 
Eu ML. Trägheitskräfte _ -agyDE-Zahl 
BE ST Schwerekräfte 
&* 

Dynamische Ähnlichkeit erfordert, daß Modell und Großausführung die 
gleiche Froude-Zahl und die gleiche Reynoldszahl aufweisen. Daraus 
folgt aber bei qy=1 (Modellfluid = Naturfluid), daß l, = 1, ist. Im 
Modellbau wird deshalb entweder Re, = Re, oder Fr, = Fr, erfüllt. 
Welche dieser Bedingungen eingehalten wird, hängt von der für jede 
Strömung individuell zu beantwortenden Frage ab, welche der drei 
Kräftearten (Trägheits-, Schere- und Zähigkeitskräfte) gegenüber 
den anderen vernachlässigbar ist. Die Froudesche Ähnlichkeit ist 
nur zu beachten bei überwiegenden Trägheits- und Schwerekräften und 
bei freien Oberflächen. Ein Beispiel dafür sind die Strömungsvor- 
gänge an einem Schiff. 

Für den Modellbau von Windrädern ist die Reynoldssche Ähnlichkeit 
maßgebend, weil es sich hier um eine Strömung handelt mit überwie- 
genden Trägheits- und Zähigkeitskräften. Die Reynoldssche Ähnlich- 
keitsbedingung lautet: 

Wurln wel 
Be HN . Gl. 34 
8 Ym YN 


RETERE 


NETTE REIT ERW DS BIER TESTER ERTENGE 


ER TGT 
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Bei wieder gleichem Fluid im Modell und in der Natur (qy=1) bedeu- 
tet das: 
Q. 


EL 
V 47 


Das heißt aber, daß die Geschwindigkeit am Modell im Windkanal 
nicht realisierbar hohe Werte annehmen müßte. Man muß sich daher 
aus praktischen Gründen mit, gegenüber der Natur, leicht erhöhten 
Geschwindigkeiten zufrieden geben. 


7. Umrechnung von Windkanalergebnissen 


ame will em Gum (mm: (Ems SEE MED siEHED AMHME SE: Mile dmmdnie sum -umiiiils SEE <mii sid AED mE si sid -mmmmuls mid MAN ME -MEME> u EHRE <mM -E EEEED UAEE ED mm 
mil EEE sin mie viällel: Gm suis Mumie SHE AMEEM MEN FE dED cöih sumiD MEN HMM sMES ED GNS AED SEM HM MEER sm: sh: smumui> dumm EEE  amEEEBS AH siMEND sEmmMER SEEN min 


EEE MUB En iD MED mEEED (ind TER ME sAdladls dm cms dm am mis wu: ai mini) su VIA WED Caihl mM MEER: MMEMEE mMED MER ER PER is EEE: SAME MED mE HER si EM sn sp +pEMEER -</MR MEE + AED mia Ad hal «chi A: ih mul cm: MED MED ME MED =MEEEB AED AMEEEE (EEEEM AED «mm 
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7.1 Motivation 


Aussagen über die Eigenschaften eines Windrades werden mit dimen- 
sionslosen Kennzahlen gemacht. Sie erlauben, mit begrenzter Ge- 
nauigkeit, Windräder unterschiedlicher Größe (siehe Abschnitt 8) 
und sogar Windräder unterschiedlicher Konstruktion (siehe Abschnitt 
3.5) zu vergleichen. Von besonderem Interesse ist die dimensionslo- 
se Leistung, der Leistungsbeiwert cp 


Die Windradleistung P wird auf die theoretische Gesamtwindleistung 


P) bezogen. 
Po=1/2*p*A, vr? Gl. 35 


Weiterhin interessiert die dimensionslose Drehzahl A,» 


A Beh. 
ct vı bl. 6) 





Die Umfangsgeschwindigkeit an der Flügelspitze ist 


uU], =ZeneneR bl. 9 
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Dies gilt auch für den Wagner-Rotor. Dabei ist R die Strecke von 
der Flügelspitze senkrecht auf die Drehachse. Der Leistungsbeiwert 
c, und die Schnellaufzahl A, werden in einem cnA-Kennfeld aufge- 
tragen. Beispiele dafür zeigt Abb. 11. 

Ar und cn berechnen sich aus Werten, die in freier Strömung leicht 
meßbar sind: Windradleistung, Projektionsfläche, Drehzahl, Dichte 
und Windgeschwindigkeit. 


Betrachtet man Abb. 8, so sieht man, daß sich bei der Abbremsung 
der Luftströmung am Windrad der Querschnitt der Stromröhre vergrö- 
ßert (Kontinuitätsgleichung). Die Vergrößerung der Stromröhre ist 
abhängig vom Abminderungsfaktor, also vom speziellen Verhalten des 
Rotors. Im Windkanal ist der Strömung jedoch an jeder Stelle ein 
Querschnitt vorgegeben. Daraus erkennt man, daß dynamische Ähnlich- 
keit (Windkanal - Natur; vergleiche 6.) nicht vorliegt, weil es 
sich um prinzipiell verschiedene Strömungen handelt. Wir wollen da- 
her versuchen, durch andere Überlegungen sinnvolle und vergleichba- 
re Kennzahlen zu erhalten. 


In Windkanal ist bei konstantem Kanalquerschnitt A, und inkompres- 


sibler Strömung auch v = const. D.h., daß man inm es, Kanal unge- 
fähr die Geschwindigkeit mißt, die auch in der Rotorebene vorliegt: 
also >E Die Geschwindigkeit v, aus Gl. 6 und Gl. 35 tritt nirgends 
auf. Eine fiktive Geschwindigkeit Vır > v, muß also erst gefunden 
werden, um den Leistungsbeiwert und die Schnellaufzahl zu berech- 


nen. 


Leistungsbeiwert c 


2 


cp 


leistungsbeiwert 





=] 
Fr 


ii 


Rotorleistungskennfeld der großen Windkraftanlage GROWIAN /2/ 


Abb. 





11 


9 = 





.  Schneraufer 





— 


FETTE: 
FFIR 

zuiEaweers 
WIRG/ZEREENAV GEHEN NG 
N 






Savonius -Rotor 


au | Leistungskennfelder 
nach Hütter /7/ für 
verschiedene Rotoran 


20 
EIRER ROSE: A 


Beispisle zum cCp-A -Kennfeld 
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42 Berechnung der fiktiven Geschwindigkeit Vır 


Die Leistungsbilanz (Bernoullische Gleichung) für die freie Strö- 


mung lautet: (vergleiche Abb. 8) 


Me 2.21 _ 
Fr z(vı v3 IE Feverluste 


Die Leistungsbilanz am Windrad für die Kanalströmung ist: (ver- 
gleiche Abb. 12) 


_M a pen 
Fk Pi Pa’r FRlerlusts 
Dabei ist (PL =P3)5 die Druckdifferenz, die durch den Rotor abgebaut 
wurde. 


Nun wird aber auch ohne den Rotor schon ein Druckverlust (Pı=P3), 
durch die Kanalwände, das Gestell und die Meßaufbauten verursacht, 
Dieser Druckverlust kann im Windkanal ohne den Rotor gemessen wer- 
den. Der Druckverlust muß bei der gleichen Windgeschwindigkeit ge- 
messen werden, bei der nachher auch der Rotor getestet werden soll, 


da (PL-Pz)y von der Windgeschwindigkeit im Kanal abhängig ist. 


(pı-pa)g*lpı"p3)-(pı-P3),, 
Gl. 36 
Ap„=Ap-Ap,, 


Die Versuchsdruckdifferenz wird im Testbetrieb des Rotors gemes- 


sen. 


Wir wollen folgende Forderungen aufstellen für unseren Vergleich 
freie Strömung - Kanalströmung: 


Forderung 1) Die Windräder sind baugleich 

Forderung 2) In der Rotorebene herrscht die gleiche 
Windgeschwindigkeit 

Forderung 3) Die Windräder haben die gleiche Drehzahl 

Forderung 4) Dem Rotor wird die gleiche Leistung zugeführt. 


We 


Wir leiten also ab: 

Die Forderungen 1,2 und 3 zusammen bedeuten, daß in der Rotorebene 
die gleichen Strömungsverhältnisse vorliegen. Es ist z.B. 
Aok "Aor (siehe auch 8.2). Die Forderung 4 wird umgeformt: 


Du SFR serR mit 61. 10 
VIF "IF olvıgtVag) P*V2r 


644- 403 Vie®Vvaett/2lvieV32) 


APR 
VIEVZEZEG VI. 


und nach Forderung 2: 


APR | 
Vie ’va.tI> Gl. 37 
IF 2K 2°p "Va 
K 
Nun haben wir aber die Schwierigkeit, daß v gar nicht exakt be- 


2K 
kannte ist, da im Windkanal durch das Windrad ein Teil der Luft zu 


den Kanalwänden abgelenkt wird (Abb. 12) und die Wind- 
geschwindigkeit in Rotorebene nicht gemessen werden kann. 


Es ist V <vV und 
2«“VYIK V2Kaußen. ’IxK 


Wir definieren: 
V=Az/A], Gl. 38 


Yist die Versperrung des Windkanals. Es ist bekannt für V>|l, d.h. 
wenn das Windrad den Kanalquerschnitt völlig ausfüllt: 


V2K”Vik 


mit Gl. 37 berechnet werden: 


ap 
FOR. 
IK Zepevig 


Es Kann damıt Vır 


SIc-Y 


Abb. 


Vak außen ? VaxK 





12 Der Wagner-Rotor im Windkanal 
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2.90: & 


für V>0, wenn also auch im Windrad Verhältnisse wie in der freien 
Strömung vorliegen, gilt: 


Yır“VYıR 


Die gesuchte Anströmgeschwindigkeit Vır soll nun durch eine Gerade 
zwischen diesen beiden Punkten angenähert werden (Abb. 13). 


V4F 


VeelV>1) 


Vae(V > 0) 
o A V 


Abb. 15 Linearisierung der fiktiven 
windgeschwindiskeit Vın 


Vır®ViplV>0)+{vi.[Ve1)-v12lV>O)}eV 


APR | 
Zen Gl. 39 





=, = 


Diese Geschwindigkeit Vir ist nun die gesuchte fiktive Geschwindig- 
keit v.; mit der die gemessenen Größen Windradleistung P und Dreh- 
zahl n dimensionslos gemacht werden können. 

Es muß noch angemerkt werden, daß bei der Umrechnung von einer 
konstanten Geschwindigkeitsverteilung über der Windradfläche ausge- 
gangen wurde und somit Störeinflüsse aus der Grenzschicht der Wind- 
kanalwände oder Verwirbelungen der Strömung am Meßaufbau nur als 
Druckverlust Ap, berücksichtigt wurden, nicht aber deren Einfluß 
auf die Anströmung der Flügel. 


7.3 Berechnung eines fiktiven Abminderungsfaktors Sr 


Um Anschluß an die geläufigen Windradformeln für die freie Strömung 
zu finden, wird ein fiktiver Abminderungsfaktor benötigt. Es ist 
jedoch zu beachten, daß dieser aus AP, berechnete fiktive 
Abminderungsfaktor einen Mittelwert über die Windradfläche dar- 
stellt und keine Bedeutung für ein einzelnes Flügelelement hat. 

Wir berechnen Sr aus Forderung 4 : 


1/2(v1p2-vig?)=Sap, 


Fa) 
V3e"YVip“Zapg 


ZAP, 


&p"V3e/Vie® rear Gl. 40 


- 32 - 


7.4 Berechnung eines fiktiven Druckes Pr 


Zur Berechnung der Reynoldszahlen Re, und Ren (Gl. 52 und Gl. 53) 
sollen zu jedem Testlauf die Größen Pı und pP, vor bzw. hinter dem 
Rotor und die Temperatur t notiert werden. In Gl. 68 (Anhang) ist 


dann für p P5x zu setzen mit 


Pı*P3-Ap,, 
Pe 2 Gl. 41 





Pr ist der fiktive statische Druck am getesteten Windrad in freier 
Strömung. 


7,5 Einfluß der Windkanalversperrung V auf den Leistungsbeiwert Co 


Der korrigierte CpnWert berechnet sich nach Gl. 1, 35 und der Kor- 
rekturgleichung 39 zu 


P 
“Peff 172°p Az Vie? 


Ohne Korrektur hätte man 


pP 
SPT7/2*p*Az*"Vv] ” 


erhalten. Der Einfluß der Windkanalversperrung wird durch den Quo- 


tienten Ap 
EEE: 3_ R 3 
Sp/Eperr"Vip"/Vig zen 
bl. 42a 


ausgedrückt. Es soll jetzt die Druckdifferenz AP,» die nur schlecht 


abzuschätzen ist, durch ee eliminiert werden: 


APR 


1 -—— Gl. 40 


sp” PVir 
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Hier wırd Vir nach Gl. 39 eingesetzt und aufgelöst: 


2 
(1-E_%)V 2 
F Ap, + l1-E.2)V-2}apotpl1-E.2)vy,2= 0 
Tepevncz R F R F K 


Die Lösung dieser in AP, quadratischen Gleichung ist: 


Zepevg.s ! 
TEE EEE ) Gl. 43 


Sinnvolle Werte erhält man mit einem negativen Vorzeichen vor der 
Wurzel. 
Das so berechnete Ap, kann jetzt in die Definitionsgleichung des 


Quotienten ep/e eingesetzt werden. Nach dem Kürzen von v,, er- 


Peff 
hält man: 


1 ET RE g 
Co/co.pgpt1t TT-E22JV [2-01-EL°)V-2V1- 01-8, ) V )}> 


Gl. 42b 


Aus Gl. 42b und Gl. 43 entnehmen wir, daß eine höhere Windgeschwin- 
digkeit zwar zu einem höheren Druckabfall führt, daß aber die Wind- 
geschwindigkeit, bei einem Windrad mit einem charakteristischen Ab- 
ninderungsfaktor er keinen Einfluß auf den Quotienten Cp/Cpaf£ 
hat. 


Die Auswirkungen der Windkanalversperrung zeigt Abb. 14. Die Be- 
rechnung des Quotienten nf Cpugr dient hier nur der Veranschauli- 
chung des Einflusses der Windkanalversperrung in Abhängigkeit von 
einer dimensionslosen Kennzahl. Bei gegebenen Windkanalmeßwerten 
erhält man diesen Quotienten natürlich schneller durch das Einset- 


zen von Ap,, vg und V in die Definitionsgleichung (Gl. 42a). 


Abb. 


14 


2,15 


2,01 
1,81 
1,43 


1,13 





Cp/Cparr als Punktion von V und g nach Gl. 42b 
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7.6 Formeln zum Versuchsbetrieb und zur Auswertung 


Zulässige Machzahl: (vergl. auch 8.2) 


Um während des Versuches Klarheit über die Machzahl zu gewinnen, 


kann man sich folgender Formel bedienen: 
| 


_ Nr VEReEeer= Eo yereerE 
M_=1/a-/(2nnR)tvr, <0,4 


Diese Gleichung erhält man aus Gl. 48 und v = vv. 


Zulässige Druckdifferenz Ap,.: 


Nachdem AP, und Ap gemessen wurden, kann Apz nach Gl. 36 berechnet 
werden. AP wird u.a. benötigt, um den fiktiven Abminderungsfaktor 
Er . Eu zu berechnen. 

E ist nur im Intervall (0;1) definiert. Mit Gl.43 kann AP, in Ab- 
hängigkeit von$ „, Yıx und V berechnet werden. Insbesondere erhält 
man mit &r = 0 und €r = 1 das theoretisch zulässige Intervall für 


R° 
2ep*vig® I 
Ep": Ap„[l&. = 0) — (2-V-2Y/Y1-V ) 
E6r1: ApnlEc=1)=0 


Abb. 15 zeigt Werte Ap,(5,=0) bei 9= S(p, und t=15"c) =1,226kg/m”. 
Diese Werte dürften theoretisch nicht überschritten werden. Die 


max. Versuchsdruckdifferenz ist 
Ap au (VıR) FAplvıg) +Apylvıg) (EL. 36) 
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20 367 336 311 290 273 
15 207 189 174 165 154 
10 92 84 78 12 68 
5 23 21 19 18 17 


Abb, 15 (Pı-P3 R bei =0 in N/m* als theorstisch 


maximal zulässigser Differsnzdruck am Rotor, 


(Yı in m/s} 
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Abb. 16 Theoretisch maximal zulässiger Differenz- 
äruck am Rotor, mit V als Paraneter 
aufgetragen über wWindzsschwindigsksit in 

"indkanal, Diagramm zu Abb. 15. 
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Größen zur Berechnung von Co und A: 


, erhalten wir aus Gl. 1, Gl. 35, P=M#2%*“*n und den Erkenntnissen 
aus Abschnitt 7. 


„4rMn Gl. 44 
p P-Azvı6” 


Mist das gemessene Moment, n die ebenfalls gemessene Drehzahl. Die 
Dichte der trockenen Luft wird nach Gl. 657 berechnet. A, ist die 
Projektionsfläche des Windrades senkrecht zur Anströwrichtung; beim 
Wagner-Rotor eine Ellipse. Es wird dabei vorausgesetzt, daß der Ro- 
tor richtig zum Wind angestellt ist, also die Achse in der Wind- 
richtung liegt. In Abb. 17 werden die Strecken und Winkel au Wag- 
ner-Rotor definiert. 





@): Umärshungswinksl 


Abb. 17 Strecken und Yinksl am Waszner-Rotor 
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Es ist 


R=a 
a3-L-sin«k 1. as 
b=L+-sink* cost ' 


Az=nab=-nL?*-sin«k- cosı 


ist die fiktive Anströmgeschwindigkeit nach Gl. 39 


Apr 
IEVIKTZeRVIg 


YıF 


.Y 


zu deren Berechnung man die gemessene Windkanalgeschwindigkeit Yıg 
® nach Gl. 67, die Windkanalversperrung V nach Gl. 38 und AP, 
benötigt. APR wiederum berechnet sich nach Gl. 36 aus dem Druck- 
verlust und der Versuchsdruckdifferenz. 
Die Blattspitzen-Schnellaufzahl A, wird aus Gl. 6, Gl. 9 und den 
Erkenntnissen aus Abschnitt 7 berechnet, 


N _ZrnR 
t Vier bl. 46 





Wir benötigen wieder die gemessene Drehzahl, den Radius aus Gl. 45 
und die Windgeschwindigkeit Yıp aus Gl. 39, 

Mit den nun bekannten Daten kann - jetzt noch einmal exakter - die 
Einhaltung der Machzahl nach Gi1.48 (Abschnitt 8.2) kontrolliert 
werden. 
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d$, Umrechnung dimensionsloser Kenngrößen 


en a I SI I ee oe 2 2 5 Io 2 5 20 85 | 


maßstäblich gieicher Windräder 


a a I a I I 4 5 4 3 7: 


8,1 Motivation 


Obwohl die dimensionslosen Kenngrößen nicht an einen bestimmten 
Rotordurchmesser gebunden sind, kann man sie, wenn eine gewisse Ge- 
nauigkeit gefordert ist, nicht von einer Baugröße auf eine andere 
übertragen. Der Grund dafür ist, daß die Strömungen i.a,. nicht 
ähnlich sind. Windräder, die die Reynoldsschen Ähnlichkeitsbedin- 
gungen erfüllen (vergl. Abschnitt 6), besitzen jedoch genau die 
gleichen dimensionslosen Kenngrößen. Es soll in diesem Abschnitt 
gezeigt werden, wie man Kenngrößen von Windrädern mit verschiedenen 
Reynoldszahlen näherungsweise ineinander umrechnet. Die Umrechnung 
basiert dabei auf dem Verhältnis der Reynoldszahlen im Modell und 
in der Natur, 
Es werden für die Unrechnung folgende Forderungen aufgestellt: 
Forderung 1} Die Windräder sind, bis auf den Maßstab, baugleich. 
Forderung 2} Die Anströwrichtung der Flügelprofile stimat 
abschnittweise überein. 
Forderung 2 kann & priori natürlich nicht erfüllt sein. Es werden 
Gleichungen hergeleitet, die dies sicherstellen. 


Einweis: Wird mit nach Abschnitt 7 aufbereiteten Windkanal- 
ergebnissen hier weitergerechnet, so ist der Index "F"” in diesem 


Abschnitt der Index "M", So ist z.B. 5, = $ 
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8.2 Die Machzahl 


Es ist nicht so sehr das Verhältnis Ma ‚/Ma, entscheidend, sondern 
sicherzustellen, daß man bei der Umrechnung von Windraddaten nicht 
Machzahlen größer als 0,4 erhält. Für Ma < 0,4 kann man die Strö- 
aung als inkompressibel betrachten, während dies für Ma > 0,4 nicht 
mehr der Fall ist. Ein Windrad in konpressibler Strömung verhält 
sich ganz anders, als eines in inkompressibler Strömung. Es ist da- 
her nicht sinnvoll zwischen diesen Zuständen umzurechnen. Da Wind» 
räder so ausgelegt werden, daß sich inkompressible Strömungen erge- 
ben und die Formeln zur Windenergie auch nur für diese Strömung 
gelten, wird für die Umrechnung die Bedingung nach Gl. 48 gefor- 
dert. 


_w 8147 
N. ao 
| 
M_=1/arvıva, ZT <0,4 1. 48 


Schallgeschwindigkeit a = [U #R&T 

a = 343 n/s für Luft bei T = 293 K 
Gaskonstante : R= 0,287 KJ/Kg K 
Isentropenexponent 9 = 1,4 


Da Ma < 0,4 an jeder Stelle des Flügels gefordert ist, muß in 
6l, 48 mit der Schnellaufzahi an der Flügelspitze Ar gerechnet 
werden. 


8.3 Schnellaufzahlen A und Abminderungsfaktoren En und Sy 


Nach Forderung 2 muß die Anströmrichtung der Profile gleich sein: 
du "Ay: Wegen Forderung 1 ist auch ß „ = ß,. Damit muß auch 
fu Py sein, weil nach Gl. 19a f =d+ß ist. Mit Gl. 18: 
tan -1/RAyo gilt letztlich: 
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Wir können Forderung 2 also wahlweise durch Gl. 49 ersetzen. 
für jeden Punkt (r,©) gilt: 





„ENT „gyag2.1/v1r17X 72.8) 
t c zZ 1/%y 1/ 1+1 Ag T1#E) bl, 29 


* 


Es ist wegen Forderung 1 


r 
M 
9] D, E, rn .. Gl. 30 und Zu" Zu 


Setzen wir Gl. 29 für Modell und Natur ins Verhältnis, so erhalten 
wir: 

C 1- 1+ 

an 5n '*®m 


NE TE TFEN 


Wir definieren 


Jrc,n/Sam Gl. 50 


nd for : 1 rn 1 - 
und formen um 301 Em)Üi+E,) (1 En li+E,) 
JEntIEn SEmen Fam en enem 
En iäpententenem  ItEmtlEn 


I-JEI-E I HEm 


ENT ITI-E HE Gl, 51 
ITI-E I +E 
Mittels En bzw. E, können jetzt Ay und A, aus Gl. 28 berechnet 


werden. Man sieht, daß gleiche Windradkonstruktion und gleiche An- 
strömung zu gleicher effektiver Schnellaufzahl A, führen; die 
Schnellaufzahlen A können aber aufgrund von unterschiedlichem Pro- 
filverhalten voneinander abweichen. Werden zwei Schnellaufzahlen 
Ay und Ax mit gleicher effektiver Schnellaufzahl Ag über Gl. 238 
berechnet, so ist die gleiche Anströmung der Flügel sichergestellt 
und danit Forderung 2 erfüllt. 
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8.4 Die Reynoldszahi 


Die Reynoldszahl ist nach Gl. 33 definiert als 





Aus den Versuchsdaten ermitteln wir die Reynoldszahl mit der fol- 
genden Gleichung, die sich aus Gl. 33 und Gl. 16 ergibt: 


Rand an R 61. 52 
eM vmlon: tm) 





Mit Vu nach Gl. 68. 

Einen gemittelten Wert der Reynoldszahl erhalten wir, wenn wir eine 
mittlere Schnellaufzahl, die an einem mittleren Radius (vergl. 
Abschnitt 83.5} gilt und eine mittlere Flügeltiefe eu einsetzen. 
Setzen wir Gl. 52 wieder für Modell und Natur ins Verhältnis, so 
erhalten wir, mit dem Geschwindigkeits- und Zähigkeitsmaßstab mit 
Gl. 17 und Forderung 2, das Verhältnis der Reynoldszahlen, das hier 


nach der gesuchten Reynoldszahl Re, aufgelöst ist. 


N 
Q, 1#E, 


=E nie (Gehege Gl. 53 
eN Rom 9,9, IHEm 


R 


8.5 Korrekturfaktoren für den Leistungsbeiwert ep 


Nachdem man durch obige Formeln in die Lage versetzt wurde, z.B. 
von einer vorgegebenen Schnellaufzahl Ay auf Ay umzurechnen, muß 
jetzt gefragt werden, wie der Leistungsbeiwert Spy Ay? aus einem 
„A -Kennfeld zu korrigieren ist, um ein Spn(Ay) zu erhalten. 

Der Leistungsbeiwert setzt sich aus mehreren Einzelwirkungsgraden 
zusammen: Dem idealen Leistungsbeiwert En; dem Profilreibungs-Ver- 
lustfaktor 7 und dem Blattzahlverlustfaktor 7 . Während der 
Blattzahlverlustfaktor für das ganze Windrad gilt, haben der ideale 
Leistungsbeiwert und der Profilreibungs-Verlustfaktor nur Gültig- 
keit für ein bestimmtes Flügelelement, Da die Vorgänge aber an je- 
der Stelle des Flügels bei beiden Windrädern ("M", "N") im gleichen 
Maßstab zueinander stehen, macht man nur einen kleinen Fehler, wenn 
man hier auch den idealen Leistungsbeiwert und den Profilreibungs- 
Verlustfaktor global betrachtet. Eine Verbesserung der Ergebnisse 
wird erzielt, wenn nicht mit den Größen an der Flügelspitze 
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gerechnet wird, sondern mit Größen, die sich auf einen mittleren 
Radius beziehen. So kann man z.B. den "Radius der halben Fläche" 


wählen: 


\ 
Ar )=1/2"Az{R) 


ır_2=1/2n7R? 
m 


rn"1/2/2-R-0,707 R 


korrekturfaktoren: 


Die Korrekturfaktoren werden berechnet, indem die Formeln zur Be- 
rechnung der jeweiligen Wirkungsgrade für Modell und Natur ins Ver- 


hältnis gesetzt werden. 


Der Korrekturfaktor für den idealen Leistungsbeiwert: 


„EPin pin An “01+2,) lo, -1) 
ae erraremeen 
cPi Cp;Mm An (rs )loy 1) 


der Korrekturfaktor für den Profilreibungs-Verlustfaktor: 
. Ren. (e nAo)liten Ag) AN, 1+Eı 
np” Npy "Te mAolklte, lo) An THE, 


Der Korrekturfaktor für den Blattzahlverlustfaktor: 


(1-1,39/(2/TP,N2))° 


NZN_ 


nz. NzMm (1-1, 39/ (zT, 2) 2 


Der gesuchte Leistungsbeiwert ist dann: 





Spy pm'Kopi Kup’ Knz 


Gl, 54 
Gl. 55 
Gl. 56° 
G1.57/ 


8 


AU 


.5 Alzorithmus zur Umrechnung dimsnsionsloser 








Kennsesrößen maßstäblich zleicher Windräder 


[Algorithmus zur Umrechnung dimensionsloser Kanngrößen | 


art ‘ 


= 


Spy! 


Profildaten: c„la,R ). ela,R,) 


Eingabe 
er Ay 9, 


fr = ,] 1 
Emm tur 1/2/2R) {'} 


= = n i 
Bm” Bm Aytr 1/2/R) {') 


nein 


ja 
Amen 


Aue, =1/2/2° 1 ch) 
MN 
Anfangswerte - 


mim Arm) 
Rum Ramt tr Am) aus Gl. 52 


A Ram’ 9, 
2 97:9, 


ara ntimal aus Profilbuch 


„Sant®-Ron) 
am'@ Ram 
Eytänlisey) aus Gl. 51 


J 


c 


it 0 CU 









Ron" Ranlay Ay: 3 Remm’ Em’ En) 


aus Gl. 53 


C (a.,R ) 

jean _ EN zus Profilbuch 
c la,R ) 
am eM 


Enäntdrtp) aus Cl. 5 





{!) 


ja 


nein 


ld) aus Gl. 27 


aus Gl. 28 
Ant ko, E,) 


4 u 
Am" 17275 im 2A 


[*} 





nein 


aralAgByn"dyn) aus Gl. 19b 










ja 
a,,g-a, l>azul. 
nein 
eyrentaR.n) 
aus Profilbuch 
Em"emla,R m) 









mais Em) 
' 





aus Gl. 26 






q 






non Ans En 





Krk p An Am Ener Em ans om) aus Gl. 55 






BEER. En. N 
Ei re a a BE 
K 









Som" Cpmirm) aus Versuch 








Spn"Spn!SpmKopi Kup-Knz) aus Gl. 57? 


Ausgabe Spy An? -ayrEnr Sans en dog"/2rro 





Für den Faktor 1/2»’/2 kann auch eine andere 
Größe gewählt werden. Hier ist entsprechend 


Kapitel B.5 ein Wert gewählt worden. Die ge- 
wählte Größe darf natürlich während des ganzen 


Berechnungsvorzanges nicht geäindert werden. 
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8.7 Anwendung des Algorithmus an einer Beispielaufgabe 


Nachdem die in den Abschnitten 8.2 bis 8.5 hergeleiteten Gleichun- 
gen in Abschnitt 8.6 in einem Algorithmus zusammengefaßt wurden, 
soll hier einmal an einer kleinen Beispielaufgabe ein Eindruck von 
der Größenordnung der Werte, insbesondere der Korrekturfaktoren, 
gegeben werden. 


Aufgabe: Im Windkanal wurden Kennfelder Say A,m? und Le 
aufgenommen. Es handelt sich um ein Windrad mit einem Profil NACA 
0015. (Siehe dazu Abb. 18 und 19). Der Flügelradius beträgt 
R,= 700mm, der Blattwinkel ist 0° und die Flügeltiefe beträgt 
t 50mm. 

Gefragt ist nun nach dem Verhältnis Son/ Spy bei einer Schnellauf- 
zahl Ay der Großausführung des Windrades (R,|=60m), wenn die Wind- 
räder bei gleicher Windgeschwindigkeit (v,=15m/s) und im gleichen 
Fluid betrieben werden. 


Bei Anwendung des Algorithmus (man vergleiche die Rechnung hier nit 
diesem) erhält man folgende Zwischenergebnisse und Ergebnisse: 


Anfangswerte: 


Atn Arm” 
An Aut 
EuAy) sei 0,4 
Re, =2,5*10 
q,”0,012 
u 
rt 
Re =2,1*10 
A=d, 
selıı 
£ =0,334 
n 


7 
) 


Rechnung: 


= * 
Re, 2,0 10 


Aog”’»#3 

Ay »15 

A =7,67° (dieser Wert ist genügend nahe am Wert Be 
um den Berechnungsvorgang schon nach einem 
Durchlauf abbrechen zu lassen) 


E_=40. En =70 


M J 


Ergebnisse: 


Das Beispiel sollte zeigen, 


ergebnisse in geeigneter Form gespeichert, 


lassen. 
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daß das Verfahren schnell konvergiert. 
Es ist in der Regel nicht mehr als ein Durchlauf nötig. Dies hängt 
natürlich von der geforderten Genauigkeit ab. 


Sind die Versuchs- 
bietet es sich an, 
ßere Datenmengen nach dem Algorithmus vom Rechner verarbeiten zu 


NACA 0019 


Test Conditions 


Reynolds Number 
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Abb. 18 Auszug aus einem Profilhandbuch /13/ 
für NACA 0015. 4,3-10%& R &3,3-10°. 


NACA 00159 


NACA VDT 
Smooth 


Test Conditions 


Tunnel: 
Date: 
Surface: 


Test: 
- Turbulence: 2.0% 


Reynolds Number 
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bb. 19 Auszur aus einem Profilhandbuch /13/ 


für NacA 0015.  6,6-10°& R,< 3,3-10°, 


Angle of Attack - "u 


Angle of Attack -a 
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9. Berechnung der Windradleistung und der 


EEE (EEE EEE mi: mi miiikie <miiiiie silmmmi: +mEEEED (AED mmmmis summml> spp ci: uimmlp (mis midi AED EHE =emiß> +miiilb (mir <A: miele SUEEElEn MEER m op sum -mpmuis mul <mmiiib ib im summum amd =EM mil: 
MEER ED -EEEED id SEA wife iii cmbiig ame ampmuifı cmmmi> ap mie SEE cmiiih mi siEN AMEED MUEEp mm -iillid: dm (ip Mining ammAMEN Amp TEEN GEMED am cm mil SED AAMED sammm cm» cm will "Pen 


Die Windradleistung berechnet sich aus Gl. |: 


P=c.*Fo 


Ist der Leistungsbeiwert bekannt, muß nur noch die Gesamtleistung 
Po gefunden werden, um die Windradleistung zu berechnen. Bei kon- 
stanter Anströmgeschwindigkeit über der Höhe durch das Windrad ist 


Po aus Gl. 35 zu berechnen. 


P, =1/2*.p-Az-vı° Gl. 35 


Da die Luftbewegung über der Erde jedoch eine Grenzschicht ausbil- 
det (Abb. 20), ist die Geschwindigkeit v, nicht konstant über der 
Höhe; sie nimmt mit der Höhe zu und kann nach der empirischen For- 


mel 


ER 61.58 


berechnet werden. Der Exponent d wird aus Messungen der Wind- 
geschwindigkeit gewonnen. Er hängt von der Bodenrauhigkeit ab. Gl. 
58 gilt nur für eine homogene Hindernisdichte auf der Erdoberflä- 
che. Im Bereich der Innenstädte gibt das Potenzgesetz nur das Wind- 
profil oberhalb der Verdrängungsdicke wieder, die etwas kleiner als 
die Höhe der Hindernisse ist. Werte für d entnimmt man Abb. 20a. 
Die Windgeschwindigkeit Vo wird in der Höhe ho gemessen. Die 


Standardhöhe ist h,=10m. 


Setzt man Gl. 58 in Gl. 35 ein, so erhält man: 


Pg=1/2=-0-Svo3Ch/ha)”dda 
A2 


= 50 = 





Abb. 20 Windhöhenprofil bei GROWIAN /2/. 


Geländeform #xponent d 


glatt (Schnse, Sand, Sse) Osclı ua 01 
leicht rauh (Grasland) 0,13: au 0,2 
rauh (Wald, Vorstadt) O2 u 02T 
sehr rauh (Stadt) O2 2a WR 


Abb. 20a Werte tür den äxponenten d aus @l. 58 /9/. 


=. 51..= 


Das Integral über die Ellipse der Projektionsfläche des Windrades 
läßt sich nicht elementar lösen. Daher wird ein Näherungsverfahren 
gewählt, welches das Integral auf eine Summe zurückführt: 
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Streifen einbeilum 
h 
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Abb. 24 Zerlegungs der Projektionsfläche des Windrades 
in horizontale Streifen, zur Berechnung der 
theoretischen Gesamntwindlsistung Po° 
m 
e 4-f: /mı2 
Po=1/2°p-»LEvı(h)?s2/Y1-(i/m)<-a-b/m 
Em t 1 IRRE. 
ST 
x Y 
Es wird definiert: 
e=b/a=b/R = = 
a Moin BL, =hz/R Gl. 59,60,61 


= coST 
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Abb. 22 Größen zur Berechnuns von ®,, 


damit: 
m 


P0.=1/2°0-1/m-2a2-ervodetg -Ch,+i-b/m) ) TÜY/TSTIZmI? 


i=-m 





Pgr1/2-0+1/m-2a2-@-Lvo3lz—gle, „-Rtire/m-R)"VT-TI/m)? 
no 


«R 


2 


",/Norenz/®n, 


‚ e m | 
Poszpatenvg’(—) l "ÜLle, +48) "O/T-Ti/m)2 


3.53 2 


Mit 


K. *f FR Gl. -63 


und 
Pu/Pm”epn/Pon’Pom/Cpm ae 3 


erhalten wir: 


Pom 


P PuKopi"Knp"Knz’p, — Gl. 64 
ON 


In Pr wird die Windgeschwindigkeit in der Höhe h, eingesetzt. Es 
ist zu beachten, daß die Höhe h bis zum Mittelpunkt der Ellipse 
geht. Dieser Punkt liegt höher als die Nabe. 


Ah=Lcos«K-sint Gl. 65 


zUI/R-Ch,, +Lcos«K-sint) Gl. 66 


eh be 


Der Faktor E50 ist im Anhang in Tabellenform gegeben. 
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10. Zusammenfassung 


Zur Korrektur der im Windkanal gemessenen Windradleistung sind fol- 
gende Schritte nötig: 

- Berechnung der Windkanalversperrung V nach Gl. 38. 

- Berechnung der Dichte nach Gl. 67 

- Berechnung des Differenzdruckes nach GLl.36 

- Berechnung des Korrektur-Quotienten nach Gl.42a 


Um eine Aussage über ein Windrad gleicher Bauart, aber anderer Grö- 
ße zu machen, folgt man dem "Algorithmus zur Umrechnung dimensions- 
loser Kenngrößen". Im Fall von Windkanalmeßwerten muß der 
ee Eur EA.) zuerst nach Gl. 40 berechnet 
werden. 


Die Windradleistung erhält man aus Gl. |: P=c ,*P. ‚ wobei die Ge- 


samtwindleistung Po nach Gl. 62 berechnet wird mit dem Faktor fg 


aus den Tabellen des Anhanges. 
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Hier im Anhang werden grundlegende Formeln zur Windenergie herge- 
leitet, die in der angegebenen Literatur nicht auftauchen, aber für 
die Berechnungen in den anderen Abschnitten benötigt werden. Wei- 
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terhin finden sich hier die Tabellen zum Faktor f, 


gramm zu dessen Berechnung. 


Gl. 16: wı Wa” u]z+vj2=vj/Yi2+1 
Gl. 27: a) Ag erhält man, wenn G1l.25 in Gl. 


eingesetzt wird 


oder: 


b) als Umformung von Gl. 
Gl. 28: Gleichsetzen von Gl]. 

ar als laser 
Gl. 35: Po=1/2mv22=0/2-V-v22=p/2-Agvp3 


Berechnune der Luftdichte o 


25 


und das Pro- 


17a und 17b liefert: 
vırad+1=1/2-v l1rg)/Ter 2 


Man betrachte die Luft als ideales Gas und wende das 


ideale Gasgesetz p=p/(RT) 


Für Luft: p9=1,293 kz/m 


„Or 


T0=273,15K 


ist der Normzustand mit: 


und 


bar 


Gl. 


67 


NT 


Berechnung der kinematischen Zähigkeit der Luft 


v=u/o u ist vom Druck unabhängig 
® uPpo 
vERsTeToI SG TeToT Don mit Gl. 67 


Linearisiert man die Kennlinie für v{pg,T) in To 


erhält man: 


v(po,T)=v(pg9,To)+1/AT-Avlpg,T)-(T-To) 


. 5 EHREER 
$ 
u er +H-rH 
Bi 
- 
” Mi 
Mi 
= = 
ZBRHERRER 
See 
Ei 
01 





au 
60 


s0 


Abb. 24 Kinematischs Zähisksit für 
atmosphärischse Luft /14/. 


v(p,T)evg*po/p+3,8°10”° m2/(s-°C)-t-pp/Pp 


merke) 


Gl. 68 


v(p,T)=pgo/p+* (vo+9,8-10. 


6 


Für Luft: vo=13,3-10 ° m2/s 


sueKe) 


sEuekz 


100 


Programm zur 3erechnuns des Faktors f20 


PROGRAM PO (INPUT,UUTPUT,TAPES=INPUT»TAPE&=OUTPUT) 


REAL N«M 
DIMENSION AusSt1il»9),EI9),ENZLLOJ,EAUS (9) 


ANFANGSWERTE 


IZ=0 
READ (5+*,END=II9) M 
DATA E/le» 0.9962, 0.9848, 0.9397, 0.8660; 
0.7660» 0.6428, 0.5» 0.3420/ 
DATA EAUS/0,,5.,10.,20.,30.,40.,50.,60.,70./ 
DATA EHZ/1le., 1.25; l.5>» 1.75» 2.» 2.5» Jo» 2 5a 6./ 


WRITE(6»100)M 
FORMATLIMI///IH »10X,' DIE TABELLEN ENTHALTEN DEN WERT 1/F0 FUER! 
„' DIE BERECHNUNG DER LEISTUNG PO.'/10%X,° 1/F9 IST ABHAENGIG *!; 
VON Ds Tau (IN GRAD) UND EHZ. M: *!,F5.) 


BERECHNUNG UND AUSGABE EINE R TABELLE 


READ (5, *+.END=9I99)N 

0021=1,10 

DOI3Y=1,9 

F=0 

DOsK=-M,M 

F=F+tEHZLII) »K/MPFELJ)) Fe L3eN)=SQARTLL-IK/M) #82) 

CONTINUE 

F=F/3.1415927/M®*2/EHZ (I) #=(3*7N) 

F=1/F 

AUStI,JI=F 

CONTINUE 

CONTINUE 

IZ=IZ+»1 

IF(IZ.EQ.4.OR.IZ.EQ.7.OR.IZ.EQ.10) WRITE(6.200) 

FORMAT IMHI//7/Y) 

WRITE(5»300) Ns tEAUSC(TIJ »I=s1»9) s {ENZCI) , (AUSLI» U) »J=1l»9) Il» 10) 
FORMAT(1INO//10X >! D: 1,F5,.3//1H »17%Xs9F7.3//10$1H »10X»10F7.37)) 
GOTO |} 


WRITE(6»400) 

FORMAT (/////* PROGRAMMENDE*) 
sTop 

END 
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IE TABELLEN ENTHALTEN DEN WERT 1/Fe9 FUER DIE BERECHNUNG DER LEISTUNG PO. 








/FO IST ABHAENGIG VON D, TAU (IN GRAD) UND EHZ. M: 100.0 
: 0.100 
N 
e>_| 0.000 5.000 10.000 20.000 30.000 40.000 50.000 60.000 70.000 
1.018 1.012 1.007 1.004 
1.01l 1.008 1.005 1.002 
1.008 1.005 1.003 1.002 
1.006 1.004 1.003 1.001 
1.00% 1.003 1.002 1.001 
1.003 1.002 1.001 1.001 
1.002 1.002 1.001 1.001 
1.001 1.001 1.001 1.001 
1.001 1.001 1.001 1.000 
1.001 1,001 1.001 1.000 
): 150 
eu rL 5.000 10.000 20.000 30.000 40.000 50.000 60.000 70.000 
1.028 1.021 1.014 1.009 1.004 
1.017 1.013 1.009 1.005 1.003 
1.011 1.009 1.006 1.00% 1.002 
1.008 1.007 1.005 1.003 1.002 
1.006 1.005 1.004 1.002 1.001 
1.004 1.003 1.002 1.002 1.001 
1.003 1.002 1.002 1.001 1.001 
1.002 1.002 1.001 1.001 1.001 
1.001 1.001 1.001 1.001 1.001 
1.001 1.001 1.001 1.001 1.000 
0200 
eu T 10.000 5.000 10.000 20.000 30.000 40.000 50.000 60.000 70.000 
1.014 1.008 1.004 
1.009 1.005 1.003 
1,006 1.00% 1.002 
1,005 1.003 1,002 
1.004 1.002 1.001 
1,002 1.002 1,001 
1.002 1.001 1.001 
1.0013 1.001 1.001 
1.001 1.001 1.001 
1.001 1.001 1.000 





=.702 


en | 0.000 5.000 10.000 20.000 30.000 40.000 50.000 60.000 70.000 





ens— | 0.000 5.000 10.000 20.000 30.000 40.000 50.000 60.000 70.000 











? 0350 
SR 

eHan 10.000 20.000 30.000 40.000 50.000 60.000 70.000 

.995 2996 998 999 1.000 

.997 2.998 .999 999 1.000 

.998 999 999 1.000 1.000 

.999 999 999 1.000 1.000 

.999 999 1.000 1.000 1.000 

1.000 1.000 1.000 1.000 1.000 

3.000 1.000 1.000 1.000 1.000 1.000 

1.000 1.000 1.000 1.000 1.000 1.000 

5.000 1.000 1.000 1.000 1.000 1.000 

6.000 1.000 1.000 1.000 1.000 1.000 





- 61 


5.000 10.000 20.000 30.000 0.000 50.000 60.000 


.973 
.983 
«988 
.992 
.994 
«996 
.997 
.999 
.999 


1.000 


.977 
.986 
.990 
.993 
.995 
.997 
.998 
.999 
.999 


1.000 


.982 
.989 
.992 
.995 
.996 
.998 
.998 
.999 


1.000 
1.000 


.988 
.992 
.995 
.996 
.997 
.998 
.999 


1.000 
1.000 
1.000 


«993 
.996 
.997 
.998 
.998 
.999 


1.000 
1.000 
1.000 
1.000 
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1. Beschreibung der Konstruktion 


Aufgabe der Sechs-Komponenten-Waage (SKW) ist, die Kräf- 
te und Momente am Windrad zu messen. Da es sich um ein 
räumlichses mechanisches System handelt, müssen sechs 
Komponenten - hier sechs Kräfts - gemessen werden. Als 
Kraftaufnehmer findet die Wägezelle 28 der Firma "Hot- 
tinger Baldwin Meßtechnik GmbH" Verwendung. Dieser Kraflt- 
aufnehmer basiert auf dem Prinzip des Dopvelbiegebal- 
kens mit aufgeklebten DMS. Die Kräfte werden über Ge- 
lenkstäbe in die Kraftaufnehmer eingelsitet. Durch den 
Binsatz. der Gelenkstäbe ist es möglich, Hystereseeffek- 


te undPassungsspiel auszuschließen. 


Das Windrad wird an der Welle des Zälektromotors befestigt. 
as wurde ein Gleichstommotor der Firma "Winkelmann" 

G100L mit einer Nennlsistung von 3 kW ausgewählt. Die- 
ser Motor/Generator erlaubt, mit seiner sroßen Leistung, 
windräder selbst bei max. Windgeschwindiskeit im wind- 
kanal, sowie Langsamläufer (hohes Drehmoment) sicher 
abzubransen. Sin weiteres Auswahlkritserium für den Mo- 
tor war die Festigksit der Welle. 3ine Überschlagsrech- 
nung ergab, daß die Welle eine ausreichende Festigksit 
aufweist, wenn die Windradunwucht nicht zu groß wird. 


Da die Kraftaufnehnmer bei Überlast leicht zerstört wer- 
den, ist die Amplitude des Krafvaufnehmerwegs durch 
einstellbare Anschläge bezrenzt worden. Die Haltsrung 
dieser Anschläge mußte besonders gsro3 dimensioniert 
werden, da sich die max. zulässigen Auslenkungen im 
Bereich von wenigen Zehntsl Millimetern bewegen. 


Der Versuchsstand kann die unterschiedlichstsen ‚indräa- 
der aufnehmen: Horizontalachsrotor, Wagner-Rotor, Darri- 
eus-Rotor. Dies wird durch sine variable Achsneizungs 
erreicht. Dazu bsfestigt man die SKY auf einem Gestell 
(Abb. 1), das im wesentlichen aus einem I-Träger besteht. 
Die unter der 5KXW windradseitig angebrachten Winkel 


werden an der obersten Bohrung des Trägers angeschraubt. 
Jeweils links und rechts verbindet ein Placheisen die 
hinteren kleinen Winkel an der SKW mit den entsprechen- 
den Bohrungen im I-Träger. Die Bohrungen sind so ange- 
bracht, daß über einen weiten Bereich eins abgsestufte 
Achswinkelverstellung mit einer Schrittweite von 5° 
möglich ist. Bei senkrechter Achsstellung fahren die 
hinteren Winksl in den I-Träger ein und werden dort 
direkt (also ohne das Placheisen) verschraubt. Der I- 
Träger wird am Drehtisch verschraubt. Durch Schwenken 
des Drehtisches ist auch eine Schräganströmungs des Wind- 
rades möglich. 


SKW wer 


r——. .——- 


Abb. 1 Skizze der SKW an Gestell mit Längenbemaßung 


2, Überschlägise Berechnuns 


Der Versuchsstand soll so Konstruisrt werden, daß der 
a-Motor voll belastet werden Kann, ohne Rücksicht auf 
die Konstruktion zu nehmen. 


2.1 Berschnung gerundetsr Maximalwerte 


Drehmoment M,: 


N/ 4 /min 750 1000 1500 2000 2500 3000 
P/yır 0,48. 9,75 1,5 250 2,5 340 
M+/Ym 6,11 7,16 9,99 9,55 3499 9,55 
Abb, 2 Tabelle zum Motor G100L: Drehzahl, Leistung 
Drehnoment 
M,=-=.10 Nm aus Abb. 2 


5 
Schubkraft auf die Windradfläche Fu: 


Bei sinem Diffesrenzdruck am Rotor, der einem Abmindse-— 
rungsfaktor & von 0,2 in freier Strömung entspricht, 
erhält man: 


"is 


P/v, 


v. = 1/2.v,.(1+8) Gl. 12 
Für v, wird die max. Windgeschwindigkeit in Windkanal 
eingesetzt, 


1/2.25m/s-1,2 = 15 on/s 


> 


F 


a 3000wW/15m/s = 200 N 


Kippmoment M,: 


Bei unsymmetrischen Rotoren, wie dem Wasner-Rotor nit 
Achs- und Konuswinkel von je 45°, tritt ein Kippmonent 
auf, welches das wWindrad nach hinten kippen will. Mit 
seinem ungefähren effektiven Hebelarm von 0,65 mal der 
Flüzellänge L ist 


M P -0,65-L 


bb” 


200N.0,65-0,7m = 91 Nm 


Unwuchtkräfte "u 


Annahme: Unwuchtmasse m = 0,015 kg bezogen auf einen 
Radius, der gleich der Flügellängse L ist. 


F e n° 


5 m. (.2°Tr) 


-L 


0,015kg+ (2.7) °-50%.(1/5)°-0,7m 


1036 N 

Achtung: Die Unwucht stellt die größte Gefahr da. Bei 
höheren Drehzahlen dürfen nur ausgewuchtete Rotoren 
benutzt werden. 


2.2 Berechnung der Gelankstäbe und der Aufnehnmer 


ee ER uhren il: ummuile Ep “eu (AUER Gm: Gem (mul (mom mem ui GE EEE MER GEmEEp min zum erumiie wind sm NEED SED AMMEp immmb veiniiis GENE Gmmma mem midi: WAR HAAN GUEE Suuimp mm 





Abb... 3 Skizze der Gelenkstabanordnung 


Stab: 93 
Der Stab 5 nimmt das Drehmoment auf, 


F, = My/r = 10Nm/0,205n = 50 N 


Stab 6: 
Der Stab 5 nimmt im wesentlichen dis Schubkraft RP. auf. 


Bei einer Vergrößerung des Achswinkels T nimmt Stab © 
auch Gewichtskräfte auf. Mit der zesamten Gewichtskraft 


F. = 800 N und dem Achswinkel T= 45° 


Rz = F+F,"sinT = 900 N 


Stäbe 1,2,3,4: 


Diese Stäbe nehmen im wesentlichen die Gewichtskraft 
N A an 17. 

F. vom Motor Fan und Rotor F„„ auf und das Aippmoment 

M,« 

FF 


G 


Betrachten wir zunächst den Horizontalachsrotor: 


Rotor 
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Anpb;- 4 Kräfte an der SKW beim Horizontalachsrotor 


Pa 2 = 1/0220, 1° (FoyFgp) = = 100 N 
P2 4 = -Fr,9 + Fan + Fan = 900 N 

Fr = PB 22 = FR, = - TON 

5 = 73 92 = R, = 640 N 


Beim ‚lagner-Rotor mit Achs- und Konuswinkel von 45° z.B. 
ist zusätzlich noch das Kippmonent M, zu berücksichti- 
sen. Jedoch verursacht die Verlagerung des Ansatzpunktss 
der Kraft Fan ein Moment in entgesengesetztsr Richtung. 
ir sehen also, daß die Abschätzung am NJorizontalachs- 
rotor in erster Näherung allzeneingültige Aussagen lie- 
fert. Durch die Unwuchtkräfts wird die Rotorkraftv en 


periodisch verszsrößert und verkleinert, 


5 = Fort Fu 


Berücksichtigt nan dies in einer Rechnung wis oben, SO 
erhält man: 


Fr} = F, zwischen -570 N und 140 N 


pn zwischen 490 N und 1060 N 


Hz) 
IN 
1 


a 


Kraftaufnehner: 


Die Kraftaufnehmer sind entsprschend den errechneten 
Kräften auszuwählen: 


Stab 5 5 kg evt. 10 kg 
Stab 6 20 kg bis 100 kg je nach Achslagzge 


Stäbs 1 und 2 100 kg 
Stäbe 5 und 4 100 kg 





Gelsnkstab: 





Abb, 5 Gelenkstab 


= F/A mit A= Deb = 12-0, 7mm“ = 8 ,Amm“ 


1000N/8,4mm® = 120 N/mm“ 


ll 


sine Aaussase uuer die Steifigkeit des Systems wird durch 


die Berschnung der Torsionseigenfrequenz gemacht. 
v=-O,um 


EEE N 





Abb... 6 Modell zum Schwingunsssystem I-Träger - Rotor- 
Motor 


n = 1/(2’T)« G-I,.(Lemer“) 
I ’ 
I, = 1/3-3 b’.s 


3s ist: m = 80 kg, G = 80000 N/mm“, ein I-Träger 
DIN 1025 T2 (IPB-Reine) 120 gewählt worden. 


I, = 1/3. (2-°120°11?+98-6,5”)mm* = 115000mn* 


damit ist n = 3,9 1/s, was als ausrsichend angese- 
hen werden Kann. 


2.4 Federkonstante von Kraftaufnehmer und Sicherung 


F 
- | LEHE 
ÜBER 
Ar f 


Abb. 7 Maße der Sicherung am Kraftaufnehner 


/(e ) 


Dir = COn; 25 ; rc 
sicherung Biege "Zug Biege "Zug 


18400 N/mm 


Die Federkonstante der Sicherung ist eine Reihenschal- 
tung aus der Federkonstanten eines Biegebalkans und 

eines Zugstabes. Die Gesamtfederkonstante von Kraftauf- 
nehmer und Sicheruns ergibt sich aus seiner Parallelschal- 
tuns. 


e =. +Ca. 
ges Aufn "Sicheruns 


Für die Hottingeraufnehmer 23 ergibt sich danit: 


F/N s/mn Sy ufn/N/am C .ag/N/um 
100 059 399 18719 
200 0,29 690 19076 
500 8,27 1852 20238 
1000 0,31 3226 21612 


Für die Grenzlast bzw. für den Grenzweg sind 1530 - 150% 
vom Nennwert angegeben. Der Anschlag soll so einsestellt 


werden, daß er bsi 100» wirksam wird. Bei Überlast wer- 
den der Kraftaufnshmer und die Sicherung gemeinsam vaer- 


fornt. Dabsi ist ein Weg s, =0,53-0,3mm zulässig. Dies 


+ 
entspricht einer Überlast P, = gas’ d+ 
F|R ist zwischen 1685 N und 1945 N. 


P_ ist abhängig vom Aufnehmer und beträst bei einem 


Grenzwes von 130% knapp 20008. 


3. Berechnuns des Lochabstandes x(®) 


Der Abstand der Bohrungen für die Laschenbefestigung 
am Gestell zur SKW für die verschiedenen Achswinkel- 
einstellungen vom Drehpunkt bsrechnet sich nach 


x(d) = L, -sinT +1 1,°-(1, .cos”T y< 


Vergleichs dazu Abb. 1. Mit dieser Formel erhäit man 
die folgende Tabelle. 


U_______ LE )/mu_ 

O 82,0 

10 1103 

’ = 

15 126,7 I, = 42 mm 
20 143,9 l, = 164 mn 
25 161,7 

30 11959 

3 197,1 
40 214,0 

45 230,1 

50 245,0 

2 250,1 

60 219,8 

65 281,3 

70 290,1 

«9 2I41 42 
80 202,0 
85 305,0 

90 306,0 bzw. direkt Verschraubung der SKW 


bsi x = 1, = 142mn 


4. Berschnuns der Steifigkeits- und der Kalibriermatrix 


Aus den ärgebnissen der Kraftmessung an den Kraftauf- 
nehmern kann man zurückschließen auf dia Kräfte und 
Momente z.B. an der HMotorwelle. Diese Umwandlung der 
Meßergebnisse erfolgt mit Hilfe der Steifigkeitsmatrix 


S. zs gilt 


wenn F der Vektor der Kräfte und Momente z.B. an der 
Motorwelle und K der Vektor der gemessenen Kräfte an 
den Krafttaufnehmern ist. Die Steifigkeitsmatrix erhält 


man, wenn bei bekannter Geometrie der SKW (Abb. 8) Kräf- 
tse- und Momentengleichgewichte aufgestellt werden. Die 
Kräfte sollen dann als positiv angenommen werden, wenn 
die Gelenkstäbe auf Zug beansprucht sind. 


M; 





Abb. 8 Geometrie und Koordinatensystem an der SKW 


F K 
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I 
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40: 2 


Mit a = 0,105 m (bis zum windradseitigen 
b= 0,317 m ände der Motorwelle) 
cC = 0,205 m 
d = 0,150 m 
ae =0,110ım 


Um die Meßeinrichtung zu kalibrieren, geht man den ent- 
gegengesetzten Weg zu oben. 

ne 

s“ 
genannt. Zur Kalibrierung wird ein bekanntss Lastkollek- 


ist die zu 5S inverse Matrix, auch Kalibriermatrix 


tiv auf die Waage aufgsbracht; die Aufnehmerkräfte wer- 
den berschnset und mit den gemessenen Werten verglichen. 


-0,340 0,350 -0,340 -1,652 0 1,652 
0,340 0,350 -0,340 -1,652 -0,16  -1,652 

71 _ |-1,047 -0,350 1,047 1,652 -4,718 1,652 

- 1,047 -0,350 1,047 1,652 0 -1,652 
O O 0 0 -4,8738 0 


OÖ -1 OÖ O O OÖ 


E-Motor 
“Winkelmann GI100L 
x 


| 


+ 
Bohrungen für alle Pußstifte 


erst bei Endmontcge ausführen 
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